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Abstract 

By including the contribution of the thermal cloud to the Lagrangian of the 
condensate of a Bose gas, we extend the time-dependent variational method 
at zero temperature to study temperature-dependent low collective excitation 
modes. A Gaussian trial wave function of the condensate and a static distri- 
bution density of the thermal cloud are used, and analytical expressions for 
temperature-dependent excitation frequencies obtained. Theoretical results 

are compared with measurements in the JILA and MIT experiments. 
PACS numbers: 03.75.Fi,05.30.Jp,32.80.Pj 
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With the development of techniques to trap and cool atoms, Bose-Einstein condensation 
(BEC) has been observed directly in dilute atomic vapors [0^01 • The new experimental 
achievements have stimulated great interest in the theoretical study of inhomogeneous Bose 
gases. The collective excitations of trapped Bose gases, which describe their dynamics and 
transport properties, have been a focus of recent experimental and theoretical studies. Low- 
lying collective excitations over a range of temperatures @-|7j and higher-lying modes || 
have been measured in recent experiments. 

In the regime with no detectable non-condensate fraction, theoretical predictions based 
on a mean-field description of weakly interacting bosons at zero temperature are in excellent 
agreement with experimental data |9|-|TT||. At zero temperature, a Bose-Einstein conden- 
sate is described by the nonlinear Schrodinger equation (NLSE), i.e. the Gross-Pitaevskii 



equation 12 



dib h 2 - , rT , . , Airh 2 a , , , 9 , 

where m is the atomic mass, U (r) the trapping potential, and a the 5-wave scattering length. 
Characterizing the evolution of the macroscopic wave function ip of the condensate, this 
nonlinear equation provides a fundamental means for studying the ground state, dynamical 
behavior and collective excitations of the condensate. The validity of the Gross-Pitaevskii 
equation in describing the evolution of the condensate wave function has been verified by 
the quantitative comparison of theoretical predictions with experimental observations. 

The Gross-Pitaevskii equation flU) has been solved by means of a time-dependent varia- 
tional technique to obtain the low energy excitations of a trapped Bose gas for both positive 
and negative scattering lengths [Q. By taking a proper trial function of a fixed shape with 



some free time-dependent parameters, a variational technique based on Ritz's optimization 
procedure leads to a set of Newton-like ordinary differential equations of the second order 
for these parameters which characterize the solution to the NLSE. Avoiding expensive nu- 
merical simulations, this method affords analytical approximations which moreover provide 
a clear physical picture of the problem. 



For a dilute Bose gas trapped in a harmonic potential, with a Gaussian function taken 
as the natural trial wave function for low energy states, expressions for frequencies of the 
collective modes in a 3D anisotropic trap have been obtained ||13|| . In the large particle 



number limit, the expressions contain the spectrum derived by Stringari ]T(J based on the 
Thomas-Fermi approximation, where the interatomic and trap interactions are dominant 
and the kinetic energy term is negligible. 

The finite-temperature excitation spectrum has also been measured in the experiments 
of Refs. || and f7|. The spectrum has been studied by means of the Popov version of 



the Hartree-Fock Bogoliubov approximation [|Tj,|l^] and by a two-fluid hydrodynamic ap 



proximation |16 , |17| . However, the finite temperature behavior of the spectrum observed in 
experiments has not been fully understood. In this paper, we extend the variational method 
at zero temperature to finite temperature where the thermal cloud exists. At a finite temper- 
ature, the thermal cloud provides an additional background field exerting on the condensate 
atoms, besides the mean-field interaction among them. The motion of the condensate can 
still be described by the nonlinear Schrodinger equation with the mean-field modified to 
include the effect of the thermal cloud. The evolution equation for the condensate may be 
derived from a variational principle to minimize the action of the corresponding Lagrangian 
density C. 

Let us consider a sample of Bosons confined in a <i-dimensional harmonic potential 

U{r) = \j:m^xl (2) 
A i=i 

where r = (x±, X2, ■ ■ ■ , Xd)- As shown in the recent experiment jTj, we may neglect the mode 
coupling between the condensate and the thermal cloud. In this simplication, we assume 
the static thermal cloud density distribution to be 

p n (r) = R(T)exp[-U{r)/kT], (3) 

where R{T) is the normalization factor. The condensate moves in the combined field of the 
external trapping potential and the mean-field potential due to both the condensate and the 



static thermal cloud. The Lagrangian density corresponding to the NLSE with the thermal 
cloud effect included can be written as 

2 \ at at ) 2m m m 

where the asterisk denotes a complex conjugate. The normalization condition for the macro- 
scopic wave function ^(r) is 

r dr\ib\ 2 = N n , 



where A<o is the number of the condensate atoms. The thermal cloud density normalizes 
according to 

J drp n (r) = N n , 

where N„ is the number of atoms in the thermal cloud. The last term in comes from the 



interaction between the condensate and the thermal cloud [18 



Following Ref. [13|, we consider the Gaussian trial function 



ij(r,t) = A(t)f[exp + , (5) 

where, by setting the Gaussian distribution centered at the origin, we have ignored the irrel- 
evant motion of the center of the condensate. The variational parameters are the amplitude 
A, width Wi and /3j relating to the curvature. All these parameters are real time-dependent 
functions, characterizing the macroscopic wave function of the condensate. Inserting the 
trial wave function into the Lagrangian density (f|) and integrating over the coordinates, we 
obtain the Lagrangian 

/oo 
Cdr (6) 
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where we have used the equation of particle number conservation 

A(t) = Nol[ 
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The last term in ([?]) is the contribution of the interaction between the condensate and the 
thermal cloud. 

By minimizing the action corresponding to the Lagrangian with respect to the parame- 
ters, we can derive the ordinary equations describing the evolution of the parameters. The 
widths of the condensate satisfy the following equation: 

w"(t) + u 2 Al - f{T)] Wi {t) = — o~~~~7 — To - + 47r f°f* 2 17 —J , (8) 

lW tl m? Wi {tf m 2 Wi(t) Y v^Fwj(t) 

where, compared with the case of zero temperature, the thermal cloud contributes an extra 
term proportional to 



UJj. 



3=1 

In the derivation we have used the fact that the widths of the condensate are much smaller 
than those of the thermal cloud, i.e. mujfwi{t) 2 <C kT. If we regard the widths (wi, u>2> ■■■Wd) 
as the coordinates of a fictitious particle, Eq. flD is just Newton's equations of motion for 
the particle in the effective potential 



1 d 

z i=i 



^ (1 - f{T))w i + 



m?Wj 2 



where the iy~ 2 term corresponds to the kinetic energy of the condensate, the /(T) term to 
the contribution from the interaction between the condensate and the thermal cloud, and 
the last a-dependent term comes from the interatomic interaction of the condensate atoms. 

After finding the widths, the curvature parameters can be obtained through the 
following equation derived from the variation principle 

mw[(t) 

m} 2h Wi (ty [ } 

Once the widths of the condensate are known, the rest of the parameters can be calculated, 
and the evolution of the Gaussian atomic cloud of the condensate can then be completely 
characterized. The whole problem of solving the NLSE reduces to that of solving the system 
of ordinary differential equations (§). 
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As a special case, we first consider a two-dimension isotropic trap. For the 'breath' mode 
with a radial symmetry, Eq. (||) reduces to 

w"(t) + u, 2 (l - h(T))w(t) = + (11) 

where, for d = 2, 

f(T) = f 2 (T) = 8n 2 N n ah\u/27rkT) 2 . 

After introducing the new dimensionless constant p = l + 2N ah 2 /m 2 and new dimensionless 
variables 

t = u^Jl - f 2 (T)t, v = w/ a , 



where ao = [h/(mJl — f2(T)u)] 1 ^ 2 , Eq. ([TT|) becomes 



v"(r) + v(r) = -^- 3 . (12) 



This equation, regarding v as the coordinate, describes the classical motion of a ficticious 
particle in the effective potential 

Vcs = v 2 +p/v 2 , 

which admits an equilibrium point only when p > 0. In this case Eq. (0) has the analytical 



solution 19 



At) - I 



v o + Z2 ± N-3 cos ( 2r ) 

v \ v 0/ 



(13) 



The oscillation frequency is 2cuyl — f2(T). For a two-dimensional harmonic trap, iV n ~ T 5//2 
21^p3|, so the temperature dependency of the frequency is 



u(T) ~ sJl-CT 1 / 2 , 

where C is a constant independent of T. At zero temperature the frequency reduces to 
2cu, as is expected. For a positive scattering length, the frequency decreases with increasing 
temperature. On the contrary, for a negative scattering length, the frequency increases with 
temperature. At a negative scattering length p becomes negative when iVo > m 2 / (2\a\h 2 ) . 



This gives a critical value of the atom number beyond which no condensate would occur, 
nor is there any oscillatory solution. This critical value is independent of the temperature. 
In the case of a three-dimensional cylindrically symmetric harmonic trap 



U(r) = - mu 2 (x 2 + y 2 + \ 2 z 2 ), 



(14) 



the equations for the widths, in terms of the dimensionless constant Po — \f^/^N Q (T)a/a 
and new dimensionless variables r = out, z/j = Wi/ao with i — x,y,z and cto 

gPza. 



dr 2 
d 2 u 



y + [l-f 3 (T)]u x 
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d 2 v z 
dr 2 



+ 



A 2 -/ 3 (T) 



y-X-Z 

v 2 z v x v y 



h/muj, are 
(15) 
(16) 
(17) 



Here, introducing the notation i] = T/T c , cu = A 1 / 3 ^; and ag 
d = 3 

f 3 (T) = f(T) = V lN n (T)- 



h/muj, we have denoted for 



(§)' 


[ ?(3) 1 




[No(T c )\ 



The fs(T) terms come from the interaction between the condensate and the thermal cloud. 
The i/j~ 3 terms correspond to the kinetic energy contribution which is neglected in the 
Thomas-Fermi approximation but plays a role in stabilizing the condensate. As in the case 
of zero temperature, for a positive scattering length a > 0, the effective potential has a stable 
equilibrium point, while for a negative scattering length a < there may be no equilibrium 
points, depending on N in comparison with a critical value N c . For an attractive interatomic 
interaction, iV c sets an upper limit of the atom number, i.e. a collapse occurs when iVo > iV c . 
When Nq < N c , a metastable equilibrium point exists as one of the two equilibrium points, 
supporting the condensate. In contrast to the two-dimensional case, for a three-dimensional 
trap, the critical atom number increases with temperature. This is consistent with what is 
known in the case of an isotropic trap |20|| . 

Following Ref. we may study the collective excitation modes of the condensate by 
expanding z/j of Eqs. (p!5|)-([T7|) around the stable equilibrium point (^ , ^o, v 0z ), which satisfies 

7 
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1 Po 



A 2 - f 3 (T) 
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(18) 
(19) 



A linear analysis then leads to the following expressions for the low excitation frequencies 



a;(|ra| 
a>(|m| 



2) = 2ujJl-f 3 (T)-2P^ 



i- 



0) = V2i 



OJ 



1 - / 3 (T)) (l + A 2 (T) )-P 2 



±J[(1 - / S (T)) (1 - A 2 (T)) + P 2 , 3 ] z + 8P| 2 



(20) 



(21) 



where = P /(4^*4) and A 2 (T) = (A 2 - / 3 (T))/(1 - / 3 (T)), and we have labeled the 
modes by the azimuthal angular quantum numbers m. By noticing that / 3 (T) = at zero 
temperature, expressions ( |20| ) and ( pi]) reproduce the results obtained in Ref. [13| . 



For small Po, expressions (|20|) and (plf) are close to the bare trap results. For large values 
of P they correspond to the case of large interatomic interactions in comparison with the 
trap excitation energies. They are good for either the repulsive or attractive interatomic 
interactions, i.e. for both the positive and negative scattering lengths. In the large atom 
number limit, where the kinetic energy term is negligible, expressions ( ^0|) and (|2lD can be 
further simplified as 



u(\m\ = 2)=uyj2(l-f 3 (T)), 



u(\m\ = 0) = uo\ 



'l-/ 3 (T) 



4 + 3A 2 (T) ± J16 + 9A 4 (T) - 16A 2 (T) 



(22) 
(23) 



which, with (1 — f 3 (t))uj replaced by u and (A 2 — / 3 (t))/(l — / 3 (T)) by A 2 at zero temperature, 
correspond to the Thomas-Fermi approximation. Near the critical temperature No is small, 
so the interatomic interaction of the condensate may not be dominant over the kinetic energy. 
For attractive interacting Bose gases of a negative scattering length, no condensate exists 
without the contribution of the kinetic energy being included. In these cases, instead of 
expressions ( p2]) and (p3"|) , expressions ( p0|) and ( |21| ) should be used. 

In the MIT experiment [7|, the axial trapping frequency v z = 16.93Hz and the radial 
frequency v r = 230Hz. Without considering the interatomic interaction, the critical tem- 
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perature kT c = hu(N\/ 1.202) 1 / 3 [|IH2|] for iV = 80 x 10 6 is 1.87//K, which is close to the 
experimental value T c = 1.7/xK. The difference may come from the neglected atomic mutual 
interaction |23|,|24[]. According to the local density approximation the chemical potential is 
given by f2~5 | 



Snapoh %uj [ 15aN 

V = + -TT \ 

m 2 I a 



1 - 



(24) 



which, for very low temperatures, reduces to 



15aN 



«o 



Noticing that a = 2.75nm in the MIT experiment, we have the chemical potential /i/^B ~ 
348nk at iVo = 15 x 10 6 in agreement with the experimental value 380nK. This verifies the 
validity of the local density approximation. 

Using the experimental parameters we have calculated the oscillation frequencies. For 
the experiment on Rb vapor at JILA, the particle number is not too large, so expressions 



(j20|) and (|2ll) are used. The temperature dependent excitation spectrum from the theoretical 
calculation together with the experimental measurements is shown in Fig. 1. It is seen that 
at low temperatures the frequencies decrease with increasing temperature. A simple theory 
explains the experimental data at low temperature when the condensate occupies a signifi- 
cant proportion. In the MIT experiment for mode m = 0, the atom number is rather large, 
so we use expressions (|22| ) and (0). The comparison of the theory with experiment is shown 
in Fig. 2. The theoretical results are in good agreement with experimental observation. 

In summary, by including a mean field acting on the condensate from the static thermal 
cloud in the Lagrangian density, we have studied the collective excitation modes of BEC 
at finite temperature by means of the time-dependent variational method, and obtained 
analytic expressions for the temperature-dependent frequency shift of low excitation modes. 
We have compared the theory with experimental observations. The theoretical results are 
in agreement with the MIT experiment very well. 
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FIGURES 

FIG. 1. Temperature-dependent excitation spectrum of the JILA experiment for collective 
modes m = (triangle) and m = 2 (circle), compared with theoretical curves (dotted for m = 
and solid for m = 2). Frequencies are normalized with respect to the radial trap frequency, and 
temperature to the critical temperature for a harmonically confined ideal gas 

FIG. 2. Temperature-dependent collective excitation frequencies of the MIT experiment for 
mode m = (circle), compared with the theoretical curve. In the experiment, temperature was 
varied by adjusting the rf frequency of a magnetic trap, so here Ai/ r f is a measure of temperature. 
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